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VIBRATIONAL RELAXATION AND RADIATIVE GAIN IN 


EXPANDING FLOWS OF ANHARMONIC OSCILLATORS 
Robert L. McKenzie 
Ames Research Center 


SUMMARY 


The equations governing the vibrational relaxation of anharmonic 
oscillators in mixtures comprised of several oscillator species have been 
solved numerically for supersonic nozzle expansions. The master equations, 
describing the rate of population change in each vibrational level of each 
species considered, include terms accounting for vibration- to-translation 
(V-T) and vibrational-to-vibration (V-V) energy transfer. A method of trun- 
cating the set of equations is developed to reduce the required number of 
equations while maintaining an accurate solution for those retained. The solu- 
tions are obtained by the use of an implicit integration scheme with the 
dependent variables represented by a set of non-Boltzmann distribution param- 
eters . The resulting population distributions among the vibrational states 
are found to display no obvious, low-order, polynomial form when large devia- 
tions from a Boltzmann population distribution are predicted. Significant 
overpopulations of the upper vibrational levels are predicted for expanding 
flows . 

The associated radiative gain coefficients are also computed for the 
infrared vibration- rotation bands of CO. They are found to be strongly influ- 
enced by the effects of anharmonic energy states on the relaxation process. 
Radiative features are sensitive mainly to the magnitude of the V-V transition 
probabilities and their variation with vibrational quantum number. The radia- 
tive gain predictions in expanding flows of CO or C0-N2-Ar are found to be of 
sufficient magnitude to support high-power laser oscillation through strictly 
gasdynamic means, even in view of the uncertainties in vibrational relaxation 
rates for expanding flows . 


INTRODUCTION 


Until recently, vibrational relaxation in flowing gases was usually 
analyzed by adopting a harmonic oscillator model of the quantized vibrational 
energy states (ref. 1) . The simple Landau-Teller dependence of transition 
probabilities on vibrational quantum number, along with the prior knowledge 
that a Boltzmann population of harmonic states will persist, then allowed the 
set of rate equations to be combined into a single expression for the rate of 
change of average vibrational energy. The model was commonly restricted 
further by considering only the transfer of energy between vibrational and 
translational modes. To verify the Landau-Teller theory, measurements were 
usually made at conditions where either vibrational excitation predominated 


(shock-wave experiments) or where only small departures from equilibrium 
existed (ultrasonic experiments). For those circumstances, the harmonic 
oscillator model provided an adequate description of the experimental results. 

Recently, however, vibrational de-excitation experiments in expanding 
flows (refs. 2, 3) have yielded seemingly large discrepancies between relaxa- 
tion rates inferred from a harmonic oscillator analysis of the expansion pro- 
cess and the rates obtained from shock-wave experiments. Those results 
stimulated further examination of the oscillator model and led to new investi- 
gations of the effects caused by small anharmonic terms in the oscillator 
potential . 

Treanor, Rich, and Rehm (ref. 4) provided a recent comprehensive analysis 
of vibrational relaxation among anharmonic oscillators, and used a "zeroth- 
order n approximation to show that the small anharmonicity of common diatomic 
gases can significantly enhance the relaxation of lower vibrational levels in 
a rapid de-excitation process. Later, Bray (ref. 5) used a first-moment 
method to obtain approximate solutions to the complete set of rate equations 
and demonstrated the manner in which upper level non-Boltzmann distributions 
develop among anharmonic nitrogen-like oscillators. 

Those recent studies (refs. 4, 5) are primarily qualitative, however, and 
apply to generalized nonequilibrium situations. Their limitations are neces- 
sary because both studies rely on some form of analytical approximation for 
solving the nonlinear differential rate equations. The work described in this 
paper was therefore undertaken to explore the results of more complete numer- 
ical solutions that not only provide a basis for evaluating the previous 
analytical approximations but, more importantly, allow some additional details 
regarding the effects of practical gasdynamic conditions to be conveniently 
included. Here, the rate equations for arbitrary mixtures of monatomic and 
anharmonic diatomic gases have been numerically integrated along the axes of 
quasi-one-dimensional supersonic nozzle expansions. The corresponding non- 
equilibrium radiative properties are also obtained since most experimental 
methods for determining the nonequilibrium state of a gas rely on its radia- 
tive properties. Furthermore, the non-Boltzmann characteristics of vibrational 
energy level populations in nonequilibrium expansions suggest several 
applications to gas laser devices. 


ANALYTICAL MODEL 


Treanor, Rich and Rehm (ref. 4) have shown that the population distribu- 
tion among the lower levels of vibrationally relaxing anharmonic oscillators 
approaches the form 


NY = N° exp 
3 3 y 



( 1 ) 


2 
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where NY is the number density of species j in vibrational state v, EY 

is the energy above the ground state, T is the kinetic temperature, and a 
is a time-dependent nonequilibrium factor, independent of species or vibra- 
tional quantum number. In rapidly cooled flows, where EY/kT greatly exceeds 
unity, av can reach values nearly equal to EY/kT, leaving the small nonlinear 

anharmonic terms in eY to describe a significantly non-Boltzmann distribution 
via equation (1) . ^ 

From the above, a tempting approximation for the anharmonic model is to 
assume distributions in the form of equation (1), and reduce the set of rate 
equations to a single expression for da(t)/dt, or some other suitable param- 
eter. The solutions to follow will show, however, that in many cases of 
interest, a varies enough with v to make such an approximation unsatisfac- 
tory for other than the lowest energy levels. In this analysis, the popula- 
tion fractions of each vibrational state are treated as separate dependent 
variables. The set of coupled, nonlinear differential equations describing 
the rate of change of each variable is then numerically integrated along the 
axis of an arbitrarily defined one-dimensional flow. No restrictions are 
placed on the distribution of oscillators among vibrational energy states 
other than to conserve their total number. Vibrational transitions are 
limited to single quantum jumps, however, since multiquantum transition prob- 
abilities are generally insignificant at the temperatures of present interest 
(ref. 6) . The populations of all rotational states are assumed to remain in 
thermodynamic equilibrium with the local translational temperature of the gas, 
and vibrational energy is allowed to be transferred only through collisions. 
(Spontaneous radiative transitions, which become important at very low pres- 
sures or for very high-lying vibrational states, are omitted.) Finally, only 
supersonic expanding flows are considered. In those cases, expansion rates 
which far exceed the collisional rates of vibrational de-excitation are easily 
obtained, and changes in the average vibrational energy of all oscillators 
remain small during the expansion, contributing little to the fluid thermo- 
dynamics. The flow properties may then be conveniently uncoupled from the 
vibrational relaxation process and the local thermodynamic state of the gas 
computed using well-known perfect gas equations for a constant isentropic 
exponent . 


The Rate Equations 


The number density of oscillators, nY, of species j in vibrational 
quantum state v may be written in terms of a mole fraction, xY, relative to 
the entire gas mixture. 


xY 


n: 




X>J 


( 2 ) 


3 


L 


In general terms, the temporal rate of change of X. due to a collisional 
relaxation process may then be written 


dxY 


dt 




,v 

jk 


( 3 ) 


where qY^ is the molar production per collision of xY due to encounters 
between species j and k, and Zj k is an average collision rate given by 


Z.. =e. d ! . 

jk m jk / 


jk 


(4) 


In equation (4), p and m are the average fluid density and molecular weight, 

2 

djk is an averaged collision cross section for encounters between species j 
and k, and yj k is their reduced mass. If the total mole fraction of species 

k is denoted as X^ = ^ \ X^, including all vibrational quantum states L, 

L 

and only single quantum transitions are considered, the vibration-to- 
translation (V - T) energy exchange contributes 


tQjdv-T - V- k 



(5a) 


while the vibration-to-vibration (V - V) energy transfer contributes 


[ Q Uv-v = v(?;»)E p jk(v;t*;) 


V v-1 V L+1 

X j h ex P\ kT 


gL+ljL E v,v-1\ 

k i 


kT 




+ p 


/L, L+l\ 
jk\v+l,v/ 


xj +1 x£ - xYx^ +1 exp 


rE L +l,L v+l,v\-] 

k i 


kT 


kT 


(5b) 


The total production term is then 


' KL * fc] v -v 


(5c) 
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Equations (5) include the detailed balance relations, exemplified by 


pV-l,V _ p V,V-l 

jk jk 


exp 



and are written in terms of relative transition probabilities 


P r > s = jk 

jk p i,o 
jk 


P 


jk 




\ 


> 


/ 


( 6 ) 


r s 

In the notation above, P_.£ is the V - T probability per collision that 
species j will make the transition from vibrational quantum state r to s 


as a result of an encounter with species k. Similarly, t ^ ie V- 

probability per collision that an encounter between species j and k will 
cause an r to s transition in species j and a simultaneous £ to m 
transition in species k. 


V 


Terms similar to 


gV , v- i 

J 


in equations (5J are the energy difference 


between two adjacent states, v and v - 1, of species j. They may be evalu- 
ated by assuming a Morse oscillator model in which the vibrational energy 
above the ground state is given by 


eY = ke.v[l - E.(v + 1)] (7a) 

j j j 

Here, 0j is a characteristic oscillator temperature and ej is the small 
anharmonic coefficient (see table 1) . 


TABLE 1.- TRANSITION PROBABILITY CONSTANTS 


Species 

Jj* 

K 

£ j 

d b 

cm 4 

CO 

3122 

6 . 203x10” 3 

1 . 36x10” 1 5 

n 2 

3395 

6 . 217x10” 3 

1 .38x10” 15 

Ar 

— 

— 

6 . 75x10“ 16 
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is then 


The energy difference. 


gV , v- i 

3 


eY ,V_ 1 = k04 (1 - 2e.v) 
3 3 


(7b) 


Transition Probabilities 


The dependence of both the V - T and V - V transition probabilities on 
vibrational quantum number continues to be one of the largest uncertainties in 
existing theoretical models.. For that reason, we shall not belabor the 
accuracy of the probability ratios Pj^ used in this study but, instead, 
explore the consequences of their uncertainties on the relaxation process. As 
Bray (ref. 5) shows, a convenient form of the theory of Schwartz, Slawsky, and 
Herzfeld (ref. 7) can be adopted in which the probability ratios are given by 



The integral expressions of the adiabaticity factors fY and , ^ 

evaluated by an empirical fit of Keck and Carrier (ref. 8) which conveniently 
bridges the gap between impulsive and adiabatic energy exchange. Keck and 
Carrier suggested the form 


,v,L 


are 


F = 




( 9 ) 


and obtained the best fit for K = 2/3. In this application, K = 2/3 is used 
as a nominal value but is also retained as an uncertain parameter. In 
equation (9) , A is given by 


A 


jk 


27r 2 £.., 






( 10 ) 


where is a range parameter describing the intermolecular potential and 

has been set equal to 0.2 A in all cases here. The magnitude of the net 
vibrational energy difference for the _V - T or V - V exchange in question is 
denoted by | AE J . The variation of P with v based on equations (8) to 
(10) and the parameters for CO from table 1 is shown in figure 1. A compari- 
son of the harmonic and Morse oscillator curves clearly shows that the small 
anharmonic coefficient, ej , has a large effect on the transition probability- 
predictions . 
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Figure 1 Normalized transition probabilities for 
harmonic and anharmonic CO. 


Equation (5) also requires 
absolute values of the ground state 

probabilities, and Pjj^i* o)* 

Since, in many cases, existing theo- 
ries do not successfully reproduce 
the absolute magnitudes of the exper- 
imental results, preference is given 
here to an empirical formulation of 
the experimental values. Data are 

1,0 

available, at least for the P 

jk 

values, because they are the most 
easily obtained from measurements of 
the net relaxation time T jk* Under 

the conditions of shock wave experi- 


ments, Tjk depends principally on 

Pjk , and may be correlated in the 
form 


jk 


B jk ( n 1 / 3 

P eX P ( C jk T 


) (ID 


It is related to 



by 


pl>0 

jk 


2 a 


8tt 


T jkP d jk X k /y jk kT 


(‘ 


-V T ) 


( 12 ) 


Values of Bj k and Cj k 
references 9 and 10. 


used in this study were determined from the data of 


Experimental data for the V - V probability P 


(°,i\ 
jkVl . 0/ 


are more difficult 


to obtain. However, what little is available and useful for this work (e.g., 
refs. 11 and 12) shows that convenient closed-form expressions like those from 

reference 7 are usually suitable for predicting the temperature dependence of 

(o , l\ /o , l\ 

P j k \l,0/ but must often be adjusted in magnitude. Nominal values of Pj^\i,o/ 

are obtained here from the following: Schwartz, Slawsky, and Herzfeld (ref. 7) 

described the transition probabilities for V - V transfer between energy 
levels far from resonance as 


(°> l ) = 
jk\l , 0 / 


0.39 tt" 


2 

P 3k 

M jMk 


^ E ik^ 

E 1 E 1 

j k 


3/2 -Oiv 
°jk 6 3 

1 _ e -< 2/3 ) a jk 


(13a) 
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where 


a =3 
jk 


2 ^kV A V 2 

h 2 kT 


11/3 


AE.. 
jk 

2kT 


(13b) 


and 


AE 


jk 



In equations (13), Mj is the reduced mass of the oscillator j and Ej 
is its energy given by equations (7) for v = 1. For the near- resonant case, 
Treanor (ref. 13) reported a modified version of the theory in reference 7 for 
small energy differences which reduces to the theory of reference 7 for exact 
resonance; namely. 


where 



u h 2 kT e 

_jk 


(14a) 


Q 


jk 


8tt 


25 A jk y jk 


( A V 2 

h 2 kT 


(14b) 


As Treanor (ref. 13) pointed out, when either equations (13) or (14) are used 
away from the region in which they apply, the prediction is too small. Hence, 
the larger of the results is used in this work. 


Truncation Methods 

With the preceding equations, all the terms of equation (3) are specified 
and the set of rate equations - one for each species j and each energy state 
v - can be written. Since chemical reactions are excluded, the set of 
equations can be reduced by one for each species by including the conservation 
condition: 


X 


j 



constant 


(15) 


The number of remaining equations is still excessive, however, since most 
oscillators of interest contain about 50 vibrational states below the dissoci- 
ation energy limit. Fortunately, if our interest is only in the lower levels, 
we can note that only negligible energy transfer occurs between the upper 
levels near the dissociation limit and those near the ground state. For many 
cases of practical interest, however, experience has shown that at least 
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20 levels must be retained to preserve accuracy in the total number of 
oscillators available. For example, if all but say the first 10 levels are 
excluded, an overpopulation of levels above the tenth could demand a large 
fraction of the available oscillators. Equation (15), summed only over the 
first 10 levels would then introduce large errors. However, even cases with 
20 levels and a single species require excessive computing times, making 
multiple-species calculations impractical. To minimize the total number of 
equations while maintaining sufficient accuracy, advantage can be taken of the 
behavior of the upper levels which couple principally through the conservation 
of oscillators expressed by equation (15) . Collisional-energy- transfer 
coupling between nearby levels, through equations (3) to (5), was found to be 
negligible beyond the nearest 2 or 3 states. Hence, little is compromised in 
the accuracy of the populations of the lower levels by properly substituting a 
partition function in equation (15) to account for those levels above the last 
one desired in the numerical integration. The rate equations are then trun- 
cated to a reasonable number (say 10 or less for each species) and only the 
population of the last few quantum states may be expected to contain some 
error. The definition of an upper level partition function requires some 
knowledge of the population distributions of those states for the anharmonic 
oscillators. The results are not too sensitive to the crudeness of the esti- 
mate, however, and a satisfactory method is to define a ’’vibrational tempera- 
ture” based on the relative populations of the ground state and the last state 
included in the integration. If we denote the quantum number of the last 
vibrational state included as v, a vibrational temperature may be defined by 



A Boltzmann distribution, described by tY is then assumed for all 

J 

states between v and the dissociation-limited last state with quantum number 
v^. Equation (15) can be replaced by 







(17) 


By approximating the values of eY in the second summation with eY 
we can obtain the numerically convenient form 



X V 

JL 

x° 


-e.v/Tj " -9jV d /T^ 

V T j 

e J J - 1 


ke .v, 
3 


CIS) 
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Mole fraction, 


Numerical experiments with the use of equation (18) showed that in the 
expanding flows of interest, truncation to as few as five levels still repro- 
duced the fourth level populations within 15 percent and had only a negligible 
effect on the ground state values. Five-level solutions without the partition 
function term in equation (18) had little similarity to more exact results. 
Some large errors still remained in the population of states near v, however, 
for cases with' developing population inversions. The inversions become arti- 
ficially enhanced by the truncation, starting at three or four levels below 
v. The difficulty was caused by a restriction in the rate equations which 
forbid transitions between oscillators in the v state and the next higher 
level. Such a restriction is a natural consequence of equation (15) whose 
differentiation leads to 


vJk - ° ™ 

k v=0 

The combined vise of equations (15) and (19) is referred to here as a "full 
truncation." The addition of v t- v_+ 1 transitions violates the condition 
of equation (19) only slightly for v > 10, and reduces the errors in XY 

near v = v considerably. Species conservation is retained by continuing 
the use of equation (15) . Hence, removal of the restriction given by 
equation (19) is an essential aspect of the truncation method and is referred 



Figure 2.— Truncation effects. (Flow conditions are defined 
by the standard case described in the text. The data 
pertain to CO at A/ A* = 200.) 


to here as a "partial truncation," 
since equation (15) is retained but 
not its derivative. 

The effects of these 
truncation methods are illustrated 
in figure 2 for a typical solution. 
The 16-level, partially truncated 
model clearly reproduces a more 
accurate 20-level model with only 
negligible errors, while the fully 
truncated 16-level model displays an 
artificial inversion. The 5-level, 
partially truncated, model also 
maintains acceptable accuracy and 
was obtained in 1/30 of the computa- 
tion time required by the 20- level 
version. Although these methods 
appear successful in the applica- 
tions shown here, their validity 
should be retested for other 
situations . 


10 



METHOD OF SOLUTION 


The successful application of many numerical integration schemes is 
often aided by the selection of dependent variables that all maintain similar 
magnitudes near unity. Such variables are usually difficult to find, espe- 
cially in chemical nonequilibrium problems where species concentrations vary 
by many orders of magnitude. Here, however, advantage can be taken of the 
nature of the solution which approaches the form given by equation (1) for the 
lower levels. The parameter a has only a weak dependence on species or 
quantum state and, for typical expansion cases, ranges through positive values 
within a few orders of unity. The a parameter therefore serves as a much 
more favorable dependent variable than XY, since the latter varies with v 

and time by many orders of magnitude. From equation (1), a convenient 
transformation is 



( 20 ) 


which allows equation (3) to be recast in terms of a\ For this application, 

the one-dimensional flow axial coordinate y shall be adopted as the indepen- 
dent variable via dy = u dt, where u is the local fluid velocity. If qY 

is defined by QY = dxY/dt = ^2 ZikQ^v* equation (3) can then be rewritten 
3 3 k 3 

for v / 0 as 


da V 

1 

dy 





kT 2 d y 


v > 0 


( 21 ) 


VO V ~ 

where X^ and X^ (v = 0) are computed in terms of and Xj 
equations (18) and (20) . 


by inverting 


The set of equations (21) (in conjunction with the fluid dynamic 
variables u, p, T, and dT/dy found from isentropic flow equations in terms 
of the area ratio A(y)/A*) was integrated as a function of - y by means of an 
implicit finite-difference scheme of Lomax and Bailey (ref. 14). Briefly, for 
this application, the method requires solutions to the linearized matrix 
equation 



where [i] is a unit matrix and Ay is the step size. The vectors W n and 
are defined at the nth step in transposed form as 
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The matrix [A n ] consists of the elements Ajj = 9Fi/3Wj at step n and 

where found by the numerical differencing procedure applied in reference 14. 
This method was found to be advantageous because it allowed the successful 
integration of the set of equations used here without further specialization 
or development. For comparison, cases computed by either a conventional or 
modified (ref. 15) fourth-order Runge-Kutta scheme developed extreme 
oscillations which eventually terminated the solution. 


Solution Characteristics 


The flow properties controlling the relaxation process are defined in 
this study by the half-angle and throat height of two-dimensional wedge-shaped 
nozzles. The starting procedure for the examples to follow assumes an initial 
value of y corresponding to a station just upstream of the sonic throat and 
a set of constant aY values corresponding to a lower state vibrational 

temperature 5 percent greater than the local kinetic temperature. These 
initial conditions approximate the progress of the relaxation process from the 
reservoir and have little effect on the downstream conditions. A standard 
case has been selected with the following general properties: 



Figure 3.— Population distribution development in CO along 
the expansion axis. Flow conditions are defined by the 
standard case described in the text. 


Nozzle half-angle 
Nozzle throat height 
Reservoir temperature 
Reservoir pressure 


<j> = 15° 
h* = 0.127 cm 
T 0 = 2000° K 
p Q = 100 atm 


The nominal transition probabilities 
referred to throughout this paper 
are defined by the parameters listed 
in table 1. 


Some resulting characteristics 
of a single-species solution are 
illustrated in figures 3 and 4. As 
expected, the vibrational state 
populations begin to distort from a 
near -Boltzmann distribution (repre- 
sented by a linear plot in fig. 3) 
as the expansion progresses, and the 
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upper levels finally become 
overpopulated to the extent that a 
net inversion develops. This over- 
population is typical (refs. 4, 5) 
when V - V energy transfer pre- 
dominates among the lower levels. 

In contrast to the behavior of 

xY(v,y), the variations of the more 
^ v 

uniform set, ctj(v,y), are shown in 

figure 4. Their difference frac- 
tions from the v = 1 value have 
been computed for the first 16 
levels and arbitrarily extrapolated 


to a . 
1 

limit. 


= 0 
Since 


at the dissociation 
v 


a , 


1 


is an exponential 


argument, even the small differences 
shown for the lower levels have a 


X v 
i * 

when the constant J a 


large effect on 


Figure 4.— Solution characteristics for a 16-level CO model. 


as evidenced 
distribution 
plotted in figure 2 is compared to 
the more exact distributions. There 
appears to be no simple, low-order 
polynomial characteristic to the 
solution, although at this time the 

waviness of the aY difference curves cannot with certainty be attributed 

solely to either a solution characteristic or a computational artifact. How- 
ever, any artificially imposed straightening caused large changes in the corre- 
sponding Xj distribution, a parameter to which the integration process is 

also significantly sensitive. The precision of the numerical method is 
implied by the scatter appearing in the most sensitive data at A/A* = 2. 

Thus, in view of the straightening effects and apparent stability of the 
implicit integration scheme, the characteristics as shown in figure 4 are not 
believed to contain significant numerical errors. 


RADIATIVE ASPECTS 


The non-Boltzmann character of the population distributions obtained in 
the preceding calculations may be expected to have some unusual effects on the 
corresponding radiative properties. Those effects have recently gained 
importance in the study of vibrationally relaxing flows because, in addition 
to their influence on the analysis of radiometric measurements, they may be 
identified as the essential feature that makes gasdynamic laser operation 
possible in diatomic expansions. 


13 



Kuehn and Monson (ref. 16), among others, obtained laser power from the 
10.6 gm vibration-rotation bands of C0 2 by placing an optically resonant 
cavity across the supersonic exhaust of a two-dimensional nozzle. .Their 
results rely on a complete inversion between the lower levels of coupled 
vibrational modes in a multimode oscillator (viz., CO 2 ) . In the previous sec- 
tion, overpopulation of upper vibrational levels in a single-mode diatomic gas 
was predicted but a complete inversion was not always achieved. Patel 
(ref. 17) has shown, however, that a complete vibrational inversion is not 
necessary for radiative gain in P-branch vibration-rotation transitions, 
provided the rotational states are sufficiently underpopulated compared to the 
vibrational states. Such a condition may be characterized by a very low 
rotational "temperature" in coexistence with a very high vibrational "tempera- 
ture." Since the non-Boltzmann vibrational distributions previously predicted 
for carbon monoxide may be represented by very high vibrational temperatures 
(as determined by the relative population of two adjacent states), and the 
rapidly relaxing rotational states follow the low kinetic temperatures of the 
expansion, we may expect to achieve radiative gain within the carbon monoxide 
vibration-rotation bands, even in the absence of a net inversion. 


The Gain Coefficient 

Radiative properties for either experimental or laser applications are 
conveniently described in terms of a small-signal spectral gain coefficient 
defined by: 


k(v) = £ [in I v ) (22) 

■where z is the optical path coordinate and I v is the spectral intensity at 
frequency v. Since laser oscillation preferentially seeks the frequency of 
maximum gain (viz., the line center), and gain measurements are most often 
made with the contribution of spontaneous emission nullified, the expression 
for gain coefficient in this application refers only to induced transitions 
at the line center frequency, v 0 . In terms of the local molecular properties, 
we can then write 


/ g u M 

k(v 0 ) = N u hv 0 B uJl n uJl (v 0 )(l - (23) 

where N u and are the number densities for the upper and lower radiative 
states, g u and g^ are their degeneracies, B u ^ is an Einstein coefficient for 

induced emission, and h U £(v 0 ) is the corresponding normalized line-shape 
function, evaluated at the line center frequency, v 0 . For purely Doppler 
broadened lines, usually considered in most laser applications, Penner 
(ref. 18) shows that 


n u£ (v o) 


2 / £ n 2 

Av d J tr 


(24) 
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where the Doppler line width is 


Av 


D 



/ 


2ftn2) (?) 


(25) 


At most of the gasdynamic conditions to follow, however, collisional 
broadening also becomes significant, and combined collision and Doppler 
broadened line profiles, evaluated at the line center, are required. The 
modified line shape at the line center then becomes 


n u* (v o) 


2 f ~i n 2 

AVp / ir 


x 2 

e [1 - erf(x)] 


(26) 


where x = /£n 2 (Av c /Avp) and the collision line width for species j is 


Av 


I £. 

tt m 




(27) 


Values of average optical cross sections, which are assumed to be 

independent of the rotational state, are listed in table 2 and reference 18 
for the gases of interest here. 


TABLE 2.- CO RADIATIVE CONSTANTS 


Bio. 

cm 2 /erg 

sec 

a 6 

02xl0 5 

Sr. 

°K 



2 

778 

6 r 




9 

05x10” 3 



( k = 

CO 

5 

3x10" 1 5 

6 jk. 

cm 2 < 


n 2 

5, 

, 3x10” 1 5 


Ik - 

Ar 

1, 

, 7x10” 15 

Based on A(v 

,j:v-l 

3-1) 

= 1.14 sec” 1 

for v 

i; j 

= 7 

at v Q 

= 2169 

cm - 1 (ref. 18) . 


Equations (23) through (27) may be specialized to the vibration-rotation 
spectra of a three-dimensional rigid rotor, anharmonic oscillator for 
vibration-rotation transitions; u -*■ 1 ~ (v, J) (v - 1 , J ± 1) where v and 
J are the vibration and rotation quantum numbers. In this case, the 
Einstein coefficients, B^, for degenerate rotation-vibration states with 
the selection rule, J = ±1, are correlated by 

V ■ B 10 V (h-Vt) < 28 > 


where 
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( 0, AJ = -1 (R-branch) 

1, AJ = +1 (P-branch) 


Since the P-branch gain coefficients are always greater, only AJ = +1 
transitions will be considered further. The molecular constant, B 10 , is 
determined from values of the commonly measured Einstein coefficients for 
spontaneous emission, A at wavelength through the relation 


B 10 




2hv 3 . v ( J * n ) 
u£ v \ 2 J + 1 / 


( 29 ) 


The ratio (g u /g & ) (N^/N u ) in equation (23) is evaluated by continuing to 
assume that the rotational states are populated according to a locally equili- 
brated Boltzmann distribution defined by the temperature T r and that 
T r = T. On that basis, we denote the mole fraction of molecules in the com- 
bined rotation-vibration state, v and J, as X v >^ and define X v = ^ X v »^ 

J 

which is identical to the vibrational mole fraction, Xj, previously computed. 

If AE V >J is the difference in energy between the Jth rotational state of 
level v and the ground rotational state of level v, then 


v /-AE V > J \ 
x exp 

gj Qr ( T r) 


(30) 


where Q r (T r ) is a rotational partition function given by Q r = T r /6 r and the 
energy difference is 


AE V,J 


k0 r J(J + 1) 1 



(31) 


In the above, 0 r is a characteristic rotational temperature defined in 
Herzberg's spectroscopic notation (ref. 19) as 0 r = hcB e /k and 6 r = a e /B e . 
The ratio gu^£,/g£^u f° r fundamental P-branch transitions is then equal to 


gj* 


v- 1 , J+ 1 


g J + l X 


V, J 


k /V- 1 


/V 


exp 


(AE V ’ J - AE V_1 ’ J+1 ) 
kT~ 


(32) 


Similarly, the upper level number density in equation (23) may be written 


N u = jjj- (2J + 1) 



(33) 
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Figure 5.— Normalized CO P-branch gain coefficients for 
X v /X v_1 = 0.8. 


The formulation given by 
equations (23) to (33) is sufficient 
to compute k(v Q ) for Av = 1 P-branch 
transitions in terms of p, T, and xY 

from the results of the preceding sec- 
tion. The necessary radiation con- 
stants are listed in table 2 for CO. 

An example of the results is 
shown in figure 5 for a gain coeffi- 
cient assuming only Doppler broadening 
and normalized by all the factors 
ahead of the brackets in equation (23) 
which are independent of J or T r . 
Note that the arbitrarily selected 
population ratio, X v /X v *" 1 = 0.8, does 
not represent a net inversion; 
although as T r decreases, positive 
gains appear and strengthen. In the 
calculations to follow, gain coeffi- 
cients are presented only for the 
value of J corresponding to the 
maximum gain in a given vibrational 
transition. 


RESULTS 


The brief studies to follow exemplify the features of vibrational 
relaxation in nozzle expansions that are chiefly an effect of molecular anhar- 
monicity or strongly modified by it. Excellent discussions of the general 
thermodynamic features of vibrational ly relaxing flows may be found in 
reference 1. 


Transition Probability Uncertainties 

In the absence of suitable experimental data with adequate sensitivity to 
the anharmonic effects, we can only indicate the limitations of our numerical 
results by first evaluating their uncertainty. The uncertainties are con- 
tained mainly in the basic vibrational transition probabilities. 

Aside from a few exceptions, most vibrational relaxation measurements 
have provided data pertaining only to the probability of transition between 
the first and ground states due to V - T exchanges (e.g., refs. 9, 10). A 
few unique investigations (refs. 11, 12) have provided additional data related 
to V - V transfer between the two lowest states of dissimilar species, but, 
in general, no suitable experimental information describing the dependence of 
probabilities on energy level for either the V - T or V - V cases is known, 
and very little data is available which defines the magnitude of even the 
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ground-state V - V rates. Furthermore, all the mechanisms governing 
relaxation in a practical experimental expansion process are not well under- 
stood. For example, some evidence (ref. 20) suggests that residual hydrogen- 
bearing impurities may cause significant increases in the "effective" V - T 
relaxation rates. In addition, recent theoretical studies (ref. 21) have 
demonstrated the importance of previously omitted long-range intermolecular 
forces at the low temperatures obtained in expansions and the interaction of 
molecular rotational excitation with the vibrational behavior. The simpler 
transition probability theories used here therefore contain intrinsic 


A Nominal rates 



:?)) 


uncertainties of significant 
magnitudes when applied to an 
expanding flow, even though 
the less sensitive shock wave 
or ultrasonic situations are 
well described by them. 

We can show the 
numerical results to be mean- 
ingful, however, if we con- 
sider only the gross features 
of the V - T and V - V 
probabilities. In the frame- 
work of this analysis, those 
features are described by the 

ground-state magnitudes P*’° 
/I 0\ 

and ^ j k \o ^ l / ’ an< ^ t ^ ie P aram “ 


(a) Population distributions. 



(b) Gain coefficients. 


eter K in equation (9) , 
which controls the slope of 
the ratios plotted in fig- 
ure 1. The influence of 
these parameters on the 
relaxation process is summar- 
ized in figures 6(a) and 6(b), 
where curves A are based on 
the nominal values stated 
earlier. The superimposed 
curves B are from V - T 
probability ratios whose 
increased dependence on 
v(viz: K = 1) leads to 

values above the twenty-fifth 
level greater than 10 times 
the nominal values. No 
effect is visible on the 
lower levels shown and uncer- 
tainties of that magnitude 
will be difficult to resolve. 
However, a factor of 5 


Figure 6.— Probability uncertainty effects for CO at 
A/ A* = 200. 
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increase in P*’° (and thereby in all the V - T probabilities , uniformly) 
has a measurable effect on all levels (curves C) . 


Variations in the V - V parameters cause different behavior. A 
decrease in the probability slope (viz: K = 1, curve D) , equivalent to the 
V - T slope change above, will be observable at levels where the gain is a 
maximum. Similarly, a factor of 10 reduction in the magnitude of p j k (o ^ l ) 

significantly affects all but the first few levels (curves E) . 

The nominal values used here were compared with experimental 

values reported in reference 12, and were found to overestimate those measure- 
ments by a factor of 10. Curves E are therefore adopted as conservative 

estimates and (1/10) ^ values are used in the remaining calculations. 

Since the V - T uncertainties have a smaller effect, their nominal values 
are retained. 


Radiative Gain in CO 

The magnitudes of gain given by curves E in figure 6(b) suggest that, 
even in pure CO, amplification might be sufficient to support laser oscilla- 
tions in an optical resonator with low losses and sufficient length. (Equi- 
valent loss coefficients are typically 0.01/L or greater where L is the 
optical path length.) Oscillation would occur preferentially at the ninth to 
eleventh levels just as found in glow-discharge CO lasers (refs. 17, 22). 

Other calculations for different gasdynamic conditions indicate that the low 
gain can be increased by several orders of magnitude by enlarging the throat 
height to 5 mm, or more. Lesser increases were obtained by halving the expan- 
sion angle or doubling the reservoir pressure. Increasing the reservoir tem- 
perature did little to increase gain but caused its maximum to shift to higher 
vibrational levels. All these trends indicate that the flow conditions chosen 
here as a standard case are, by no means, optimum and appear to be restricting 
the gain due to excessive freezing of the V - V, as well as the V - T, rate 
processes . 

The radiative gain computed from a harmonic-oscillator model has also 
been included in figure 6(b). For these conditions, the anharmonicity terms 
clearly modify the distribution of gain within the vibration-rotation band. 
Except for predicting radiation from the lowest level, the harmonic-oscillator 
model is unsuitable for any radiative analysis of this type. 


Vibrational Exchange in Mixtures 

When a second vibrationally relaxing gas with different relaxation rates 
and similar energy level spacing is added, vibrational energy is exchanged 
between the two species and enhances the vibrational excitation of one at the 
expense of the other. The addition of N 2 to CO expansions has a particularly 
large influence on the behavior of CO because of two features of N 2 : (1) the 
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vibrational energy levels of N 2 are nearly coincident with those of CO, 
allowing rapid V - V energy transfer between the two; and (2) the V - T 
relaxation rates of N 2 are slower than those of CO, allowing N 2 to act as an 
energy reserve for CO in expansions while making vibrational energy available 


2 . 0 - 



to CO through V - V collisions. 

The vibrational "pumping" of CO by 
N 2 can cause the equivalent CO 
vibrational temperature to exceed 
even its initial reservoir tempera- 
ture. This behavior is illustrated 
in figure 7 for vibrational tempera- 
tures described by the first and 
ground state populations (see 
eq. (16)) of each species in several 
CO-N 2 admixtures. 

Others (refs. 4, 23) have shown 
that if V - V transfer between the 
two species is rapid compared to the 
V - T rate, the vibrational states 
of the two are related by 


Figure 7.— First level (v = 1) vibrational temperatures in 
CO-N 2 mixtures. ( v — V probabilities are 1/10 the 
nominal values.) 


9 C0 

fTv Co\e 



(TV, 


N 2 qe 


6 C0 - 6 n 2 
- 


(34) 


without regard to the anharmonicity of either species. In equation (34), 
(Ty„_) is the first-level vibrational temperature for a quasi-equilibrium 

condition maintained under the requirements just described. In expansions of 
anharmonic oscillators, equation (34) is satisfied initially but becomes less 
representative as the expansion progresses. A comparison of Ty^ obtained 
from the numerical calculations and from equation (34) using numerical values 
of Ty. and T is shown in figure 8. The disagreement that develops between 

the two is caused by freezing of the V - V process. The later freezing in 
mixtures of greater CO content is due to the larger V - V probabilities for 
CO-CO encounters than for C0-N 2 collisions. 
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Figure 8.— Deviations from quasi-equilibrium V — V 
coupling in CO-N 2 mixtures (1/10 Py _ y). 



I 


Gam Enhancement in Mixtures 

The enhancement of CO vibrational excitation caused by V - V exchanges 


with N 2 also leads to increased gain as 



Upper level vibrational quantum number, V 


Figure 9.— Gain coefficients in CO-N 2 mixtures at 
A/A* = 200 (1/10 Py_ v ). 



Figure 10.- Gain coefficients in Ar-CO-N 2 mixtures for 

^N 2 /^CO = 19 (1/10 Py- v). 


figure 9 illustrates. The CO popula- 
tion distribution, in the presence of 
N 2 , becomes more uniform (higher 
equivalent TV C q) and more CO mole- 
cules are placed in the upper levels, 
even though the total number density 
of CO is reduced. Both factors con- 
tribute to increased gain. Figure 9 
also indicates that gain is not gen- 
erally optimized by any single ratio 
of C0:N 2 . The mixture giving maximum 
gain depends on the vibrational tran- 
sition of interest as well as the 
gasdynamic conditions. The general 
conclusion may be made, however, that 
the addition of N 2 over a wide range 
of mixtures raises the CO gain coeffi- 
cients well above the equivalent opti- 
cal loss coefficients conveniently 
obtained in most laser cavities. 

Gain may be enhanced by other 
than V - V pumping as well. Fig- 
ure 5 indicates that gain may also be 
increased by reducing the rotational 
temperature. This must be done with- 
out affecting the degree of vibra- 
tional excitation, however (i.e., 
without reducing the reservoir tem- 
perature) . Since rotational relaxa- 
tion requires only a few collisions, 
the rotational energy or temperature 
in most expansions usually remains in 
thermal equilibrium with the local 
translational temperature. At a 
given location in the expansion and 
from given reservoir conditions, both 
temperatures may then be reduced 
through strictly thermodynamic means, 
by adding large fractions of a mon- 
atomic gas. Argon is particularly 
suitable for this purpose because, as 
a heavy collision partner, it also 
reduces the rate of vibrational 
energy loss through V - T colli- 
sions. An example of the net effect 
is shown in figure 10, where an N 2 /C0 
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mixture ratio of 1:19 has been diluted with argon. The gain is substantially 
increased and its maximum is shifted to lower levels. 




Diatomic Gasdynamic Lasers 


5 . 3 - 

P(I6) 


k(i/ Q ) 


The applications of laser devices depend, in part, on the wavelengths of 
maximum gain or power. In a diatomic gasdynamic system, the vibrational and 
rotational level of maximum gain may be shifted downward by going to larger 

area ratios where the rotational 
temperature is lower; the gain is 
increased, and in some nozzles, the 
optical path is lengthened. This 
would give rise to laser oscillation 
at shorter wavelengths with increased 
power and offer some wavelength selec- 
tion. Some examples of these varia- 
tions are shown in figure 11 for an 
expansion that is limited to 
A/A* = 4000 by the conditions for 
condensation of CO to a liquid state. 
While the lower level gain (v = 1) has 
Figure 11.— Wavelength and gain variation along the reached a maximum determined by the 
expansion axis for a 25% CO-75% N 2 mixture decreasing gas density, higher levels 

(1/10 Py„y). are expected to reach maximum gains 

limited by other physical properties not included here, particularly gain 
depletion due to spontaneous radiative decay. In general, the operation of 
powerful, gasdynamic, diatomic laser devices appears feasible, even in view of 
the large relaxation rate uncertainties associated with expanding flow 
predictions. 1 



Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, October 23, 1970 


■^Since the completion of this theoretical study, the operation of a CO 
gasdynamic laser based on the principles described here has been demonstrated 
(ref. 24). Its qualitative behavior is in accordance with these predictions. 
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